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ABSTRACT 

Context. Whilst there is a generally accepted evolutionary scheme for the formation of low-mass stars, the analogous processes when 
moving down in mass to the brown dwarfs regime are not yet well understood. 

Aims. In this first paper we try to build the most complete and unbiased spectroscopically confirmed census of the population of 
Collinder 69, the central cluster of the Lambda Orionis star forming region, as a first step in addressing the question of how brown 
dwarfs and planetary mass objects form. 

Methods. We have studied age dependent features in optical and near-infrared spectra of candidate members to the cluster (such 
as alkali lines and accretion associated indicators). In addition, we have complemented that study with the analysis of other youth 
indicators like X-ray emission or mid-infrared excess. 

Results. We have confirmed the membership to Collinder 69 of ~90 photometric candidate members. As a byproduct we have 
determined a temperature scale for young M, very low-mass stars and brown dwarfs. We have assembled one of the most complete 
Initial Mass Functions from 0.016 to 20 Mq. And, finally, we have studied the implications of the spatial distribution of the confirmed 
members on the proposed mechanisms of brown dwarfs formation. 

Key words. Stars: formation - Star: low-mass, brown dwarfs - open clusters and associations: individual; Collinder 69 



1. Introduction 

In the current paradigm, stars are bom within molecular clouds, 
which are accumulations of gas and dust. These clouds are ini- 
tially supported against gravitational collapse by a combination 
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tational instabilities yields a collapse of the cloud material ( Shu 
|et al.|198"7l and references therein). 

The discovery of the first brown dwarfs in 1995 (very low 
mass objects characterized by the lack of stable hydrogen burn- 
ing in their interior, masses typically below 0.072Mq) led to a 
debate on the formation mechanism of this type of objects that is 
still open today. Since the typical thermal Jeans mass in molecu- 
lar cloud cores is around 1 M©, a thermally supported cloud does 
not fragment in cores of substellar masses and the formation of 
brown dwarfs cannot be directly explained as a scaled-down ver- 
sion of low mass star formation. 



While Padoan & Nordlund ( 2002 1 argued that brown dwarfs 



form via "turbulent fragmentation" (the density enhancements 
produced by the turbulence decrease the Jeans mass); Reipurth & 



Clarke ( 2001 1 suggested that they may be stellar embryos ejected 



mass to start hydrogen burning (in this model, the truncation 
of the disks is explained by dynamical interactions). Besides, 
more recently, Whitworth & Zinnecker ( 2004 1 proposed photo- 



from newborn multiple systems before they accreted enough 



evaporation of massive pre-stellar cores as the formation mech- 
anism of brown dwarfs. 

Even though a significant number of brown dwarfs with 
ages around a few Myr have been reported to harbour active 
disks (|Luhman et al.||1997 Fernandez & Comeron 2001 , Natta] 
'et al. 2004; Barrado y Navascues & Mar tin 2003; Barrado^ 
iNavascues et al. 2004a Mohanty et al. |2005) favouring the 
"in-situ" formation scenario, more homogeneous and system- 
atic studies of disk accretion, rotation and activity (and the re- 
lationship between classic indicators of activity such as Ha and 
these phenomena) in young brown dwarfs could help to confirm 
the "universality" of this mechanism or the dependence on other 
(environmental) factors. 

These kinds of studies can only be carried out setting their 
bases on robust census (well below the hydrogen burning limit) 
of star forming regions of different ages and environments. This 
work presents such a census for a very interesting star forming 
region, and, by studying the spatial distribution of its members, 
addresses the question on the mechanism of formation of brown 
dwarfs on this particular environment. 
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The Lambda Orionis Star Forming Region (LOSFR) is asso- 
ciated to the 08III star A Orionis (located at ~400 pc,|Murdin & 



Penston 1977), the head of the Orion giant. It comprises both re- 
cently formed stars up to ~24 Mq and dark clouds actively form- 
ing stars. Although its properties (morphology, distance, redden- 
ing, size) make this star forming an ideal laboratory to test star 
formation theories, until recently, it had not been very well stud- 
ied. 



At the beginning of the 80s, Duerr et al. ( 19821 carried out 
an Ha emission survey identifying three stellar clusters centered 
around the dark clouds Barnard 30 and Barnard 35. Those clus- 
ters were later confirmed from a statistical point of view by 



Gomez & Lada| ( [1998) . In particular, Collinder 69, the central 



one, is a well-defined, compact open cluster affected by rather 
low extinction Av~0.36 mag, Duerr et al.|1982) . It is quite rich, 
with one O binary star (A Ori itself), about a dozen B stars (Duot| 
|et al.|1982[ ) and a well-populated sequence of low-mass stars and 
brown dwarfs ( Barrado y Navascues et al.|2004b [2007) . 

More recently a number of photometric studies have been 
published focused on Collinder 69: Dolan & Mathieu (1999 
[2001 1 [2002| obtained optical photometry and presented a selec- 
tion of candidates with estimated masses down to ~0.3 Mq (as- 
suming an age of 5 Myr). jMorales-Calderon ( 2008| l presented a 
very complete compilation of photometry from the optical to the 
mid-IR obtained with different ground-based and space obser- 
vatories (including previous candidates by [Barrado y Navascu6s[ 
et al.|[200 4c 2007). They propose a list of candidates with es 



timated masses well below the hydrogen burning boundary if 
membership to the cluster is confirmed. This sample (together 
with the X-ray candidates, see next section) represents the start- 
ing point of this work. Our aim here is the spectroscopic char- 
acterization and confirmation of membership. When we com- 
bine our membership results with those of other spectroscopic 



2004c 



works ( [Dolan & M athieu 1999 2001 ; Barrado y Navascues et al 
Sacco et al. 2008 



tion rate of only - 



Maxted et al. 2008) we find a pollu- 
this suggests that the methods followed by 



Morales-Calderon ( 2008 1 to obtain candidate members are very 



reliable. 

This paper is organized as follows: In Section[2]we describe 
very briefly the photometric and X-ray data from which the can- 
didate selection was obtained by Morales-Calder6n[ ([2008 ) and 
Barrado et al. ( 201 l| l. In Section|3[we describe the spectroscopic 



observations of the candidates. In Section |4] we present a tem- 
perature scale derived for M young very low mass stars and 
brown dwarfs (from our spectroscopically confirmed members 
of Collinder 69). In Section |5] we describe the procedure fol- 
lowed to build our final census of members. Finally in Sections|6] 
and[7|we present the study of the spatial distribution of the con- 
firmed members (and its implications in the theories of forma- 
tion of brown dwarfs) and one of the most complete Initial Mass 
Functions constructed so far, respectively. Our conclusions are 
summarized in Section[8| 

On a forthcoming paper (Bayo et al. 2011, in prep., from now 
on Paper II), we will present properties of individual members of 
Collinder 69 (such as accretion rates, rotational velocities, etc.) 
as well as more general cluster related ones (such as disk ratios 
and disks spatial distribution). 

2. Photometric candidates in Collinder 69 

As mentioned before, the base of this work is the follow-up 



of the photometric candidates proposed in Morales-Calderon 
( 2008| l and a sample of the X-ray emitters reported in Barrado 



|et al.| ( [201 l| l. Here we provide a short description of the data 



used in those two works and in Fig [T] we show a diagram with 
the coverage of the different surveys over the contours of the 
lOOyum IRAS image of the cluster: 

- Optical: 



1. 



The CFH T1999 Survey (Cousins R and I bands): de- 
scribed in Barrado y Navascues et al. (2004c i, for cluster 
members, the faint limit is set by Rcompiete~22.75 mag at 
(R - I) = 2.5, corresponding then to Icompiete,duster~20.2 
mag. For a DUSTY 5 Myr isochrone (Chabrier et al. 
2000), and the distance and standard extinction for the 
cluster, this limit corresponds to 20 My„p„(,r. 
The Subaru2006 Survey (i' and z' bands): described in 
Morales-Calder6n[(p008]l, for a DUSTY 5Myr isochrone 



( [Chabrier et al.|2 000), and the distance and standard ex- 
tinction for the cluster, an I magnitude of 25.5 mag cor- 
responds to ~7 M.jupiiei-\ for a COND 5Myr isochrone, 
more appropriate for this range of temperatures, the com- 
pleteness limit is located at 3 M-jupner- 
Near Infrared: 

1. 2MASS (J, H, Ks bands) provides near infrared data 
down to a limiting magnitude of 7=16.8, H-lbA, and 
££=15.3 mag (-30 My„,„-,„. for a DUSTY, [Chabrier etal. 
[2000[ 5Myr isochrone). 

2. INGRID Survey (J, H, Ks bands): described in [Barrado 



y Navascues et al. ( 2007 1, the detection limit can be es- 



timated as J/,,„„f=21.1 mag and the completeness limit 
as Jcompfete=19.5 mag. For a DUSTY 5 Myr isochrone 
( [Chabrier et al.|[2000| l, and the distance and standard 
extinction for the cluster, this limit corresponds to 

10 M.j„pi,er- 

3. Omega2000 Survey (J, H, Ks bands): described in 
[Barrado y Navascues et al. ( 2007| l, the completeness lim- 
its for flie survey are: JcompZc/c=20 mag, Hco,„,,fe,c=19, 
and Kscomptoe- 18. For a DUSTY 5Myr isochrone 
( [Chabrier et al.|[2000 ^, and the distance and standard 
extinction for the cluster, this limit corresponds to 

8 Mjiipiier- 

Mid Infrared: 

1. Spitzer mid-IR imaging (3.6, 4.5, 5.8, 8.0 and 24 . yum): 
described in Barrado y Navascues et al. (2007i 



and |Morales-Calder6n ( 2008 1, the completeness lim- 
its are at [3.6]compiete=16.5 mag, [4.5]compiete=16.5, 
[5.8]compiete=14.5, and [8.0]compiete=13.75. For cluster 
members, the completeness limit at 3.6 fim (for a 5Myr 
isochrone by Baraffe et al.|1998" I corresponds to a mass 

-0.04 Mq. 

X-rays: 

1. XMM-Newton observations of Collinder 69: The obser- 
vations consisted of two fields, one to the East and an- 
other to the West of the bright star A Ori and allowed 
us to study the Weak-line TTauri population of this open 
cluster. A list of detections was compiled, cross-matched 
with our previous photometric surveys and several se- 
lection criteria were applied to derive a catalog of new 
candidate members (see details of the whole process in 
[Barrado et al.[201 l] l. This study produced a list of 66 can- 
didates (19 new) from which 44 sources have been spec- 



troscopically confirmed already by Dolan & Mathieu 
([ I999| l; [Barrado y Navascues et al][ ( |2004c[ ); [Sacco et al. 
|2008t ; paxted et al.| ( [2008| l or this work (4 of them have 
only been confirmed by us). The X-ray detected clus- 
ter sample is complete down to ~O.3M0, with some 
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detections for confirmed members with masses close to 
O.1M0. 



10.5 



9.5 




R.A. (deg) 

Fig. 1. Sketch of the areas covered by the main optical, infrared, 
and x-rays surveys in the region of CoUinder 69. A Ori itself is 
highlighted with a large star and the B population with smaller 
ones. Omega 2000 and Ingrid data together fill the area covered 
by the CFHT survey (we do not include the mosaics for the sake 
of clarity in the figure). The large red dashed boxes show the 
IRAC FoV -the the edges are the regions with coverage only 
in 3.6/5.8 or 4.5/8.0- and the area surveyed by MIPS is slightly 
larger. 



3. Spectroscopic observations and data analysis 



Regarding the spectral analysis we have used our own data 
(see Section s |3.1| and 3.2 



et al. (2008 ); Maxted et al 



and the measurements from 



Dolan 



& Mathieu' (2001 1; [Barrado y Navascues et al] ( |2004c| l; 



(12008]). 



Sacco 



In brief: Ipolan & Mathieu (2001 ) obtained R~20000 specti-a 
with a setup that allowed them to study the Li I /16707 A and 
Ha lines, and a rich array of metal lines near 6450 A for pre- 
cise radial velocity measurements to confirm members down to 
~0.3 Mq. Strong Ha emission (indicative of accretion) has also 
been detected from three candidate substellar members in the 



very deep survey for young st ars of iBarrado y Navascues et al. 
( |2004c| l. Finally, [Sacco et al.| (|2008|) and [Migted et al.| ( |2008| l 



published their results on the comparative analysis of high reso- 
lution spectra of two samples of low-mass young stars (candidate 
members to cr Ori and Collinder 69 clusters); with a resolution of 
R~16000-17000 they studied binarity, Li absorption. Ha emis- 
sion and derived rotational velocities for candidates from D olan] 
& Mathieu ( 1999 1 and Barrado y Navascues et al. (2004c| close 
to the brown dwarf domain. 



Because of the large amount of data we were dealing with (in 
particular remarkable for own our data-sets, see Tables [T] and |2] 
and more details in Bayo 2009 ), we developed an automatic pro- 
cedure to perform the line characterization. We also studied the 
effect of the spectral resolution in the measurements provided by 
this automatic procedure (see Appendix A). 

3.1. Optical spectroscopy 

During the last seven years, our group has been granted time 
at different observatories to perform spectroscopic observations 
of the previously described candidates. These observations com- 
prise a wide range of resolutions and wavelength coverages. In 
Table [T] the most relevant information of the different runs is 
displayed. 

Except where otherwise stated, we reduced all data within 
the IRAhM environment in standard way. 

LRIS, Low Resolution Imaging Spectrometer (Keck): We 

used the 1200 lines/mm grating, with a scale of 0.63 "/pixel and 
a resolution of ~2 A (measured in a NeHe comparison lamp, 
R~2650). The one arsec slit was used and a wavelength coverage 
of 6425-7692 A was achieved. During the same run we collected 
low-resolution spectra with the 400 1/mm grating and also the 
one arcsec slit. The wavelength calibration is better than 0.4 A, 
the resolution is 6.0 A around the wavelength of Ha as measured 
with a NeAr lamp (R~950), and the wavelength coverage: 6250- 
9600 A. 

MIKE, Magellan Inamori Kyocera Echelle (Magellan): 

We used the prior to 2004 Red-MIKE standard configuration: 
R2 echelle grating, scale of ~0.29"/pixe and 1.0" slit. With our 
set-up, the spectral coverage in the red arm was 4430-7250 A. 
In order to improve the final S/N, we degraded the resolution by 
rebinning the original data during the readout to 2 and 8 pixels 
in the spatial and spectral directions, respectively, achieving a 
resolution of 0.55 A (R~11250). 

B&C, Boiler & Chivens spectrograph (Magellan): We 
used the 1200 1/mm grating and the 1.0" slit; with a scale of 
0.79"/pixel a resolution of ~2.5 A as measured in a HeNeAr 
comparison lamp (R~2600); and a wavelength coverage from 
6200 to 7825 A. We also used the 300 1/mm grating (same slit as 
before, a resolution of ~800), and a wavelength coverage 5000- 
10200A; this way, the Magellan spectra have slightly worse res- 
olution and larger spectral range. 

FLAMES, Fibre Large Array Multi Element 
Spectrograph (VLT): Program 080.C-0592. We used the 
LR6 grating with a measured resolution of 0.76 A; each fiber 
has an aperture of 1.2 " on the sky and GIRAFFE has a scale 
of 0.3"/pixel in MEDUSA (the resolution is -8600 and the 
wavelength range 6438-10350A). The data reduction was 
performed using the GIRAFFE gir-BLDRS pipeline vers. 1.12, 
following the standard steps which include correction for the 
differences in the fiber transmission. The spectra processed by 
the pipeline are not corrected for sky background; therefore, 
the analysis that we present was performed after subtracting 
a sky spectrum. This background spectrum was computed as 
the median of those obtained from the fibers positioned "on 
sky". Since the nebular emission of the region of Collinder 
69 is not negligible and our sky fibers were distributed quite 



' IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc. under contract to the National Science 
Foundation. 



3 



A. Bayo et al.: Spectroscopy of Very Low Mass Stars and Brown Dwarfs in the Lambda Orionis Star Forming Region. 



homogeneously, we studied the variations of the Ha nebular 
emission with these sky fibers to have an idea of the accuracy of 
the correction achieved. We measured a mean full width at 10% 
of ~41 km/s with a standard deviation of ~3 km/s. Therefore the 
dispersion measured in different fibers translated into an added 
~7% uncertainty in our measurements. 

CAFOS, Calar Alto Faint Object Spectrograph (CAHA 
2.2m telescope): We used the R-200 grism, with scale of 
0.53"/pixel and a measured resolution of ~11 A (we used the 
1.6" slit, R~600, with a wavelength coverage of 6200-10350A). 

TWIN (CAHA 3.5m telescope): In every run we used the 
T-13 grating and the 1.2" and 1.5" slits and only the data coming 
from the red arm were processed. The pixel scale is 0.56"/pixel, 
and we measured a resolution of ~6 A (R~ 1 100 and wavelength 
coverage from 5600 to 10400 A). 

3.2. Near-infrared spectroscopy 

SOFI, Son of ISAAC (NTT): Program 078.C-0124. We used 
two low-resolution grisms with the 0.6" slit to roughly cover the 
JHK bands on a Hawaii HgCdTe 1024x1024 detector with a 
plate scale of 0.292"/pix. The blue grism covers 0.95-1.63 fim 
and the red grism the region between 1.53-2.52 jum. The cor- 
responding spectral resolutions were 930 and 980 respectively. 
The telescope was nodded along the slit between two positions 
following the usual ABBA pattern. 

In addition to the science targets we observed several "tel- 
luric standards" (AOV objects at similar airmasses) to remove 
telluric water absorption bands as described by |Vacca et al.| 
( 2003| l, and to estimate the instrumental response. A Xe arc lamp 
was used for the wavelength calibration (consistent with the OH 
airglow calibration) with an accuracy of 1 .2 A for the blue grism, 
and 2 A for the red one. 

NIRSPEC, the Near InfRared echelle SPECtrograph 
(Keck): In both campaigns we used Nirspec-3 with the 0.57" 
slit, covering the 1.143-1.375 fim wavelength region. The cor- 
responding spectral resolution is ~2000 (and the pixel scale 
0.18"/pixel). Again, classical ABBA pattern was followed and 
telluric AOV standards were observed. In this case, the data re- 
duction process was carried out using the IDL based software 
REDSPEC (which performs the standard steps in an pseudo- 
automatic manner where some interaction with the user is re- 
quired). 

IRCS, the InfraRed Camera and Spectrograph 
(Subaru): We used the "Grism HK" covering a wavelength 
range of 1.4-2.5 fim at a very low resolution of ~150 (0.3" 
slit and the 52 mas scale). The spectra of the science targets, 
spectral templates and AOV standards were obtained following 
an ABBA pattern and the reduction of the data was performed 
using IRAF in a standard manner. 



4. Temperature Scale for young M very low-mass 
stars and brown dwarfs 

4.1. Spectral Typing 

Depending on the expected nature of the sources themselves and 
the characteristics of the available spectra, we have used differ- 
ent approaches to derive spectral types for our candidates: 

Optical spectra: We considered two groups according to the 
effective temperature derived from the SED fit (see Sect. |4.2[ i: 
"warm" (Teff > 4000 K, about MO spectral type) and "cool" 
(Teff < 4000 K) sources. 



For the warmer part of the sample (a fraction of the XMM 
candidates), we compared our optical (low resolution, CAFOS) 
spectra with templates (obtained with the same configuration 
during the same campaigns) from Taurus members and field 
dwarfs within a spectral type range from GO to MO (the compar- 
ison spectra were reddened to the previously mentioned average 
low extinction. Ay ~ estimated by |Duerr et al.|1982| for Collinder 
69). We normalized both, the science spectra and the templates, 
at the same wavelength and by a simple minimization decided 
which template reproduced the science data best. We did not use 
specific temperature sensitive lines since our resolution did not 
allow us to fit anything else than the continuum shape. As an ex- 
ample, in Fig. [2] we show an example of this comparison for a K 
candidate member. 

On the other hand, the majority of our candidates belongs 
to the colder sample, in principle with Teff consistent with M 
spectral types. M-dwarfs are characterized spectroscopically by 
the presence of molecular bands of titanium and vanadium oxide 
(TiO, VO). These bands reach their maximum strength around 
M7, and then become weaker for cooler temperatures because of 
Ti and V condensation into dust grains. Regarding atomic lines, 
the Ca II triplet is much weaker than in M supergiant spectra, but 
on the other hand, the Na I and K I doublets are stronger, since 
they are gravity dependent. 

Several spectral ratios or indices (quantifying different band 
strengths) have been proposed in the literature to derive spectral 
types for these late-type stars. Some are based on the relative 



depths of the mentioned molecular bands (see for example Reid 
|et al.||T995j |Cruz & Reid 2002 and references therein); others 
are based on measuring the slope of the pseudocontinuum (for 
example PC3 and PC6 from , Martin et al. 1 1996, and |Martm et al] 
[T9991 ). 

We have used different combinations of these indices (de- 
pending on the resolution and wavelength coverage of the spec- 
tra) to classify our M-like candidates. In Fig.|2]we show an ex- 
ample of a spectral sequence obtained applying these indices to 
13 candidates observed with FLAMES. This combination of in- 
dices should provide us with a spectral classification with ~0.5 
subclass accuracy. 

Near-infrared spectra: Because of their nature, the reddest- 
coldest candidates of our sample (late M, L and even T spectral 
types according to the SED fitted temperatures), can be better 
studied in the infrared than in the optical. As in the optical case, 
we performed the spectral classification either through spectral 
indices (for the two M candidates, see right panel of Fig. [3]l or 
by comparison with templates (for the L and T candidates, see 
left panel of Fig. [3]l. 

The near-infrared spectra of M dwarfs are dominated by deep 
broad absorption bands of H2O (particularly at 1.4 and 1.85 fim). 
The fact that water vapour in the Earth's atmosphere also con- 
tributes with a substantial absorption at these wavelengths ham- 
pers the analysis of these bands to some extent. However, the 
higher temperatures in the stellar atmospheres mean that the as- 
sociated steam bands are broader than the terrestrial absorption, 
so the wings are accessible for measurement and analysis. 

Apart from these water bands, other molecular bands such as 
the CO (at 2.29 fim), FeH (at 0.99 yum ) and VO (at 1.2 fim) also 
scale with temperature (Jo nes et al.| [l994). Nevertheless, these 
changes are not as dramatic as the ones observed on the water 
bands and, therefore, we will focus on the former to derive the 
spectral types. 



Comeron et al. ( 2000 1 defined a reddening independent in- 



dex, /H2O, to measure the depth of the wings of the water band 
of late M dwarfs centered near 1.9 yum. Gomez & Persi ( 2002| l 
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Table 1. A summary of the Collinder 69 optical spectroscopic campaigns. 



Date 


Observatory/Telescope/lnstrument 


Resolution 


Wavelength 


Number of sources 


Original photometric 






AA/A 


coverage 


observed 


survey 


Nov 2-5, 2002 


Mauna Kea / Keck / LRIS 


-2650 


6425-7692A 


12 


CFHT1999 


Nov 2-5, 2002 


Mauna Kea / Keck / LRIS 


~ 950 


6250-9600A 


29 


CFHT1999 


Dec 11-14, 2002 


Las Campanas / Magellan / MIKE 


-11250 


4430-7250A 


14 


CFHT1999 


Mar 9-11,2003 


Las Campanas / Magellan / B&C 


-2600 


6200-7825A 


2 


CFHT1999 


Mar 9-11,2003 


Las Campanas / Magellan / B&C 


-800 


5000-10200A 


3 


CFHT1999 


Nov 22-25, 2005 


CAHA / 3.5m / TWIN 


-1100 


5600-10425A 


5 


CFHT1999 


Nov 20-23, 2006 


CAHA / 3.5m / TWIN 


-1100 


5700-9900A 


8 


CFHT1999 & 1" X 1° Spitzer 


Nov. 30 -Dec. 11,2007 


CAHA /2.2m / CAFOS 


-600 


6200-10350A 


37 


CFHT1999 & 1° X 1° Spitzer & 








6438-7 184A 




XMM-Newton survey 


Jan 5, 2008 


Paranal / VLT /FLAMES 


-8600 


40 


CFHT1999 & 1° X 1° Spitzer 



Table 2. A summary of the ColHnder 69 near infrared spectroscopic campaigns. 



Date 


Observatory/Telescope/Instrument 


Resolution 


Wavelength 


Number of sources 


Original photometric 






AA/A 


coverage 


observed 


survey 


Dec 22-23, 2004 


Mauna Kea / Keck / NIRSPEC 


-2000 


1. 143-1. 375;um 


4* 


CFHT1999 


Dec 9, 2005 


Mauna Kea / Keck / NIRSPEC 


-2000 


1. 143-1. 375yum 


9* 


CFHT1999 


Jan 9-11, 2007 


La Silla / NTT /SOFI 


-950 


0.950-2.500yum 


2 


CFHT1999 


Nov 10, 2008 


Mauna Kea / Subaru /IRCS 


-150 


1.400-2.500;um 


8* 


Subaru2006 



* For some of the sources the obtained S/N ratio was too low due to poor weather conditions to perfonn the analysis and therefore are not listed on the tables with 
the results from these analyses. 
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Fig. 2. Left: Spectral type determination of one XMM candidate 
member by comparison to templates obtained with the same in- 
strumental set-up (CAFOS). The "science spectra" are displayed 
in black and the templates (labelled according to the spectral 
type) in red. Right: Spectral sequence derived for candidate 
members of Collinder 69 observed with FLAMES. Note the in- 
tense Ha emission and strengthening of the typical TiO and CaH 
bands with the spectral type. 



tested this method by comparing the estimated spectral types 
with the Q index defined by [Wilking et al. ( 1999| l and obtained 
excellent agreement (see Fig.8 in their paper). This Q index is 
reddening independent too and was defined in order to character- 
ize the strength of the 1.7-2. 1 fim and > 2.4 fim water absorption 
bands. 



Because of the wavelength coverage of our observations, we 
only used the /h,o index to estimate infrared spectral types of the 
candidates: In the right panel of Fig. [3] (dotted line) we show the 
linear relation proposed by Comeron et al. ( 2000) between this 
index and the spectral type (obtained for confirmed members of 
Chamaeleon I, ~1 Myr old). We have re-calibrated this relation 
using a sample of well-known field M-dwarfs (filled black cir- 
cles with error bars illustrating the dispersion in the measure- 
ments made on ~3 objects per spectral type) obtaining a very 
similar slope (red line). With this comparison, we see how, in 
the M5.5-M9 spectral range, the /hjO index is not sensitive to the 
known age dependency on the water bands. Keeping this caveat 
of the age uncertainty in mind, this relationship allows us to es- 
timate spectral types within 1.5 sub-type. For the two objects 
where we used this method we also had optical spectra. In both 
cases, the infrared spectral type is colder than the optical one but 
within the error bars. This might be indeed a result of the age de- 
pendency of the /h,o index. Finally both sources were classified 
as diskless based on their IRAC photometry, so the possibility of 
some excesses caused by the disk affecting our classification can 
be ruled out. 

Both optical and near-infrared spectral types derived by us 
are listed in Table |6] 



4.2. SED fitting 

As can be inferred from Section[2j the photometric surveys pro- 
vided us with a very large but inhomogeneous dataset. To ana- 
lyze these data in a homogeneous and automatic manner we used 



the tool VOSA ( |Bayo e t al. 2008). 

In short, the user provides VOSA with a table with pho- 
tometry and the tool performs the following tasks in an almost 
automatic way (very little interaction is needed): it enlarges 
the wavelength coverage by looking for counterparts in differ- 
ent (VO-compliant) catalogues, it builds the SEDs for all the 
sources, it performs the SED fit using mo dels ([H auschildt et al. 
IT9991 lAllard et al.||200T] [Chabrier et aLllSOOOi ,Castelli~erar 
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Spectral Type 



Fig. 3. Left: Infra-red spectral type determination for a candidate 
member observed with Subaru. The templates used for compar- 
ison (red spectra with corresponding spectral type label) were 
also obtained with Subaru and correspond to known L and T 
dwarfs from Geballe et al.|2002 A visual comparison with the 
spectra obtained by |Lodieu et al | ( |2008[ ) of Upper Sco members 
also suggests L spectral type as a better match specially based 
on the water bands and the "peaky" structure ~L7yum. Right: 
Our own re-calibration of the reddening independent /hjO index 
from Comeron et al. (2000^. The dotted black line represents the 
linear relation found by Comeron et al. ( 2000 1, while the red line 
is the best linear fit to our data and is the relation that we have 
used to estimate spectral types of our M candidates. 



1997\ of stellar and substellar photospheres (effective temper- 



ature, surface gravity and metallicity as free parameters), it uses 
those models to obtain a multi-color bolometric correction and, 
finally, it interpolates among isochrones and evolutionary tracks 
( Baraffe et al.||1998| [20021 [Chabrier et ai:||2000| [Baraffe et al.| 
[2003 ) to derive ages and masses. 

In addition to the photometric data, the user has to provide 
VOSA with the estimated distance to each source and the ex- 
tinction affecting the line of sight. In the case of Collinder 69, 
this does not represent any extra complication, since members 
should be located approximately at the distance derived for the 
cluster (400 pc, Murd in & Pens ton[[ 1977 1 and we benefit from 
the low extinction (Av~0.36 mag, Duerr et ar]|1982 1 affecting 
this region. 

The parameters derived in this fashion are listed in Table [6 
The relation between this VO methodology and other classica 
ones can be found in [Bayo et al.[ ( [2008) 1. In particular, some of 
the possible caveats to take into account in the determination of 
Teff, Lboi and Mass via SED fit are: 

- Blue excess in the SED due to accretion. This aspect is 
discussed in Paper II: there we study in detail the veiling in the 
objects with the strongest Ha emission and demonstrate that for 
these cool objects veiling is hardly noticeable and therefore the 
determination of Teff, Lboi and Mass via SED fit is not sensitive 
to it. 

- Red excess in the SED due to the presence of a circum- 
stellar disk. The impact of circumstellar disks on our SED fits 
is discussed in [Bayo et al.| (2008), where we showed that only 
edge-on disks above a certain mass (below which, the disk is not 
massive enough to produce an effect in the photospheric part of 
the SED) introduce a bias in our parameter estimations. When 
this is the case, our methodology will underestimate both the ef- 
fective temperature and the bolometric luminosity of the central 
object. From the SED shape of the confirmed members (no in- 



frared peak higher than the photospheric one) we can infer that 
we do not have edge-on disks in our sample. In any case, to illus- 
trate this point more clearly, in Fig. [4 we plot the estimated spec- 
tral type vs effective temperature. While the latter could in prin- 
ciple be affected of the presence of the disk, the former should 
not since we are using blue features of the optical spectra to es- 
timate the spectral type. In the figure we have highlighted the 
objects that harbour disks by plotting large circles around them. 
We can see how the dispersion in effective temperature among 
the objects of the same spectral type does not correlate with the 
presence of disks. 

- Gray excess in the SED due to multiplicity. To study the ef- 
fect of possible unresolved multiple systems we have performed 
a simple exercise considering two cases: a similar mass ratio and 
an "extreme" one. In both cases we have built simulated SEDs 
by adding the fluxes of two pairs of confirmed members and 
performed the fit in the composite SED. As expected, when the 
components have similar characteristics (in this particular case 
M4 spectral type, diskless sources); the derived effective tem- 
perature does not change, but VOSA over-estimates the mass 
and Lboi of the individual components (and of course, underes- 
timates the age). On the other hand, for the most extreme case 
that we could build (class III sources but with 5000 and 2800 K 
effective temperatures), the resulting Tgff is very close to that of 
the hottest component (4700 K) and the Lboi larger than that esti- 
mated for the hottest component, but still almost within the error- 
bars. Thus, in this case again, the larger effect to keep in mind 
is that the age of the system would be underestimated. In any 
case, according to [Sacco et"ari ( |2008| l and [Maxted et al.[ ( [2008] l, 
the, -short period-, binarity fraction in Collinder 69 is as low 
as ~10% (unlike the higher ~30% reported for field M dwarfs 
by [Reid & Gizis[|1997] l; so, although for the individual targets 
we should keep these possible biases in mind; the general con- 
clusions about the cluster as a whole should not be affected by 
binarity. 

4.3. Effective temperature scale. 

Several temperature scales for M dwarfs are already available in 
the hterature ([Bessell et al.[[l998; Luhman 1999; Leggett et aL 



2000 Luhman et al.[[20()3^ However, in many cases, the sam- 
ple of objects for which this scale is derived lacks homogeneity 
and/or is based on small numbers statistics. As an example, the 
set of objects studied sometimes comprises sources with differ- 
ent ages or environments, and therefore the possible effects of 
those factors cannot be addressed. 

With this in mind, using the determinations of spectral types 
and effective temperatures from the previous two subsections, 
we have derived our own temperature scale. We are taking ad- 
vantage of the homogeneity of our determination of spectral 
types and effective temperature and of the fact that we are study- 
ing objects with roughly the same age, that were born in the 
same environment (an environment with the practical advantage 
of low extinction affecting the observational data). To derive the 
temperature scale, we started by defining a "clean" sample (see 
Table [3] and Fig. [5]l. This "clean" sample is composed of ~30 
sources fulfilling every one of the next criteria: membership con- 
firmed through lithium absorption (see next section), SED with 
no infrared excess whatsoever (to avoid problems in Teff deter- 
mination implied by edge-on disks with VOSA, see Section 4.2 
and Bay o et al.|200 8), optical spectral type derived between MO 
and M8 (see Section [4T[ ), no signs of variability or binarity and 
estimated age according to the HR diagram between 1 and 10 
Myr (see Fig.[TT). 
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Table 3. Temperature scale derived from our specttoscopic data 
of CoUinder69. 



1000 7 



2000 : 



3000 7 



4000 7 



5000 : 



Spectral Type Effective Temperature 
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4000 K 


ML5 


3750 K 


M2.5 


3600 K 


M3.5 


3500 K 


M4 


3500 K 


M4.5 


3300 K 


M5 


3200 K 


M5.5 


3260 K 


M6 


3100 K 


M6.5 


3050 K 


M7 


3000 K 


M8 


2700 K 



6000 i 



MO 



M5 

Spectral Type 
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Fig. 4. Comparison between the spectral types determine from 
spectra and the effective temperature estimated via SED fit. 
We only display the spectroscopically confirmed members here, 
highlighting the sources that harbour disks with large circumfer- 
ences. It is clear that no systematic differences in Teff estimation 
can be attributed to the presence of disks. 



The temperature scale determined in this manner is shown 
in Table [3] I n Fig. [5] we co mpare our scale with others from 
the lit erature ( |Luhnian|1999) [Bessell et al.|1998| [Luhman etal] 
2003[l. W e find a good agreement with the scale derived by 
Luhman (1999i for a sample of Taurus members (1-3 Myr, 



slightly younger than Collinder 69) but systematically higher 
temperatures than those proposed by Bessell et al. ( 1998 1 and 
Luhman et aL] ( |2003| l. With respect to the [Bessell et al.| ( |1998] l 
scale; what we see is probably just the effect of the difference 
in age from the two samples of objects for which the scales are 
derived. 



(2003 



On the other hand, the scale presented by Luhman et al. 
correspond to an update of the one derived in Luhman 
( |1999| l, but in this case for the young cluster IC348 (~2 Myr 
old). Therefore, the differences between the two scales in this 
case should not be related to an aging effect. In any case, these 
differences are not significant taking into account the dispersion 
in estimated T^ff that we see for an individual subclass (specially 
between M3.5 and M6). 

In this same line, we must note that the dispersion previ- 
ously mentioned, seems to be of the same order as the differences 
found by Luhman ( 1999[ l between dwarfs and giants (while none 
of our confirmed members can have a low luminosity class ac- 
cording to our study). This dispersion seems to be higher within 
specific spectral types and might be revealing some physical 
properties of the atmospheres of objects within this range of 
temperatures. It could be related to dust settling, although the 
temperature range seems to be too high for this phenomenon to 
happen. It could be related to metallicity, but, to the best of our 
knowledge, there are no detailed metallicity studies in this clus- 
ter In any case, our "clean" sample is not well enough sampled 
to allow us to derive any conclusion. 
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Fig. 5. Temperature scale derived for the confirmed very low 
mass M star and brown dwarf members of Collinder 69. With 
different symbols and line styles we also display scales previ- 
ously reported in the literature. 



5. Confirming membershiip 

To confirm the membership of our candidates we have used dif- 
ferent diagnostics (all of them based on alkali absorption lines) 
depending mainly on the resolution of the spectra: 



5. 1. Lithium absorption at 6708 A: 

Above ~0.065 Mq (less massive objects cannot develop the nec- 
essary temperature in their cores), lithium acts as an age scale 
because the time it takes for the core to reach 3.0x10^ K is a 



sensitive function of mass (Basri 1997 1. Very low mass stars and 
brown dwarfs down to this mass limit are fully convective, and 
therefore, once the core temperature exceeds the necessary limit, 
the entire lithium content of the star should be exhausted rapidly 
and thus be reflected in an observable change in the photospheric 
lithium abundance. Theoretical models make specific predictions 
about the time evolution of this lithium depletion boundary. For 
example, Ventura et al.| ( [T99 8 ) predict that at ages 30, 70, and 
140 Myr, the lithium depletion edge should occur at 0.17, 0.09, 
and 0.07 Mg, respectively. Other models by Chabrier & Baraffe 
( 1997\ and Burrows et al. ( 1997 1 make similar predictions of the 
variation of this lithium depletion boundary with age. Indeed, 
Bildsten et al.| ( |1997) and others have argued that the age for 
an open cluster derived in this manner should be more accurate 
than those found by any other method ( [Stauffer et al.|1998|[T999l 
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Barrado y Navascues et aL|1999 Jeffries et aL|2003 [Barrado y 
Navascues et al.||2004c[ IJeffries & Oliveira||2005t |Jeffries et aL 
2009l[ 



Because of the intrinsic weakness of the lithium line we have 
only been able to use it as a youth indicator for the sources ob- 
served with a high signal-to-noise ratio (especially in the lower 
resolution campaigns) and down to a certain spectral resolution 
(see Appendix A). Whenever these two factors were favorable, 
the lithium presence (and absence) was considered the main 
component to assess membership in Table |6] As an example. 
Fig. |6] shows two objects with clear Li I absorption in spectra 
taken with different instrumental set-ups. 



LCri043 U4.Q 



C59WI1 -2708 M3-S 



6b7G 




)690 6700 6710 6720 6720 67§670 668: 
Lombda (A) 



Fig. 6. Detail around Li I /16708 A for eight confirmed members 
with diff'erent spectral types. The two panels correspond to two 
different campaigns: Magellan/MIKE echelle spectra (only the 
order corresponding to this line), and FLAMES multi-fiber spec- 
tra. 



To illustrate the presence of Li I in absorption as a youth 
indicator. Fig. |7] shows the measured lithium equivalent width 
versus the spectral type of the (newly confirmed) members to 
Collinder 69. In the middle panel, we have included the data 
from |Dolan & Madueul ( [19991 |200T) ; [Sacco et al.| ( [20081 ) along 



with the measurements from this work. Whenever more than one 
measurement of the Li equivalent width was available we display 
the average value. For comparison, in the bottom panel of the 
same figure, we have included data corresponding to a similar 
age cluster, cr Orionis, as well as (in all the panels) an upper 
envelope of the values measured in older clusters (see caption 
of the figure). We have also plotted the theoretical equivalent 
widths from Zapatero Osorio et al.j (p002 ) considering log g - 
4.0 and \ogg = 4.5 and initial cosmic abundance (A(Li) - 3.1). 

The scatter of the Li I equivalent widths is considerable for 
all spectral types included in the figure, getting even larger in 
the M domain. This is a well-known trend for which there is not 



a clear explanation yet (see 


Barrado y Navascues et al.||2001a| 


Zapatero Osorio et al. 2002 


and more recently da Silva et al. 


2009 and references therein). The dispersion could be ascribed 



to a variability in the Li I line as a consequence of stellar activ- 
ity, diff'erent mixing processes, presence or absence of circum- 
stellar disks, binarity, or different rotation rates from star to star 
Recently, jBaraffe & Chabrierl ( |2010| l tried to explain this disper- 



sion in terms of the early accretion history that could lead to a 
dispersion in the Li abundances of young low mass stars and 
brown dwarfs. In the top panel of Fig. [7| we have plotted red 
dots on top of the sources in our sample classified as accretors 
(based on the Ha equivalent width, see Bayo et al 201 lb, here- 
after Paper II, for details), large open circles around those show- 
ing infrared excess and large open squares for those sources clas- 



sified as binaries by either Sacco et al. ( 2008 ) or Maxted et al. 



(2008 1. It turns out to be quite clear that none of these special 
sets of objects show distinct positions in the diagram. Therefore 
the cause for the scatter observed in these measurements remains 
unknown. 

Based on this scatter, [da Silva et ari ( |2009l ) warned about the 
possible problems of using Li to date individual stars; we have 
tried to quantify whether this can be a problem in our case, and 
we have found that in all the cases (16 sources) where member- 
ship was confirmed by Dolan & Mathieu (2001) or this work 
based on Li and where Maxted et al.| ( 12008 ) had studied the ra- 
dial velocity; membership was confirmed by them too. On the 
other hand, we do not have means to measure the opposite ef- 
fect, since for the only two objects that we have discarded as 
members based merely on lithium (LOriOl 1 and LOri012), there 
is no radial velocity measurement to compare with. In any case, 
for very young associations like this one, and for the estimated 
spectral types of these two objects (K8 and K9), even with scat- 
ter, lithium should be detectable (see Fig 
[2009l for more details). 

In addition to the dispersion in the lithium EWs among 
sources with different temperatures, for some objects we have 
observed some variability in this line when comparing our data 
with the values measured by jSacco et al.j ( |20081 l or Dolan &1 



2 of |da Silva et al. 



Mathieu] ( [19991 [200 l| l (see also Sec 5 on Paper II). We display 
our measurements against theirs in Fig. [8] It looks like most of 
our values are in good agreement (within the error-bars) with 
those of the literature. However, for the sources listed in Table[4j 
the differences cannot be explained in terms of uncertainty of the 
measurements. Most of these objects show features typical of ac- 
tive stars (such as X-ray emission or Ho' variability). Neuhaeuser] 



et al. ( 1998 ) monitored the young star Par 1724 finding a variable 



lithium equivalent width consistent with rotational modulation. 
In any case we only have two measurements per source, and the 
sample is too small to reach any conclusion. 

5.2. Sodium and potassium absorptions: 

Regarding the lower resolution sample of sources, other alkali 
lines (more prominent than Li I) are known to be gravity sensi- 
tive in M-type stars. This is the case for some sodium and potas- 
sium doublets: K I at 7665 & 7700 A and Na I at 8184 & 8195.5 
in the optical and K I at 1.169 & 1.177 /urn and 1.244 & 1.253 /urn 
in the near-infrared (see jMartin et al.|1996[|Schiavon et al.|1997| l. 
The surface gravity of members of Collinder 69 in this range of 
temperatures is expected to be log(g) = 3-4 ( [Baraff'e et al.j 
1998 ), whereas a typical M-type giant will have \og{g) ~ 2 and 
an older main sequence dwarf will have log(g) ~ 5. Thus, the 
equivalent width of these doublets can be used to identify (and 
discard as members) background giants and foreground dwarfs 
from our sample. 



To illustrate these differences, in Fig. 12 we show a compar- 
ison of the Na I doublet of candidate sources with field dwarfs 
of the same spectral types. 

Due to the inhomogeneous nature our data, we had to chose 
the most suitable lines to be compared with templates or mod- 
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Table 4. Properties of the sources showing hthium EW variability ( Sacco et aL 2008[ Dolan & Mathieu 1999 



2001 values compared 



to the values derived in this work, third column). Note that there is no bias concerning the resolution of the spectra (as was also 
proved with the exercise of Appendix A). The instrument code is the same as in Table|6] 



Object 


SpT 


INS 


Binarity' 


IR Class 




vsini' 


Ha Vble* 


X-ray^ 


LOri013 


M3.5 


(8) 




III 


3750 






Y 


LOri045 


M3 


(8) (6) 




III 


3500 


<17 


Y 


N 


LOri057 


M5.5 


(3) 




III 


3700 


<17.0 


N 


N 


LOri063 


M4.0 


(3) 




II 


3700 


<18.2 


Y 


N 


LOri068 


M4.5 M5.0 


(1)(3) 




III 


3700.0 


<17.0 


Y 


N 


LOri075 


M5.0M5.5 


(2) (8) (3) 


SBl (M08) 


III 


3400 


61.3!iV' 


Y 


N 


LOri088 


M5.5 


(8) 




III 


3200 


<17.0 


Y 


N 


LOri094 


M5.5 


(8)(1) 




III 


3200.0 




Y 


N 


LOril06 


M5.5 


(1) 




II 


3200 


<17.0 


Y 


N 



Maxted et al. ( 2008b (M08). 



' According to'Sac co et al.| ( |2008| (SOS) or 
2 Deriv ed with VOS A. 
^ From |Sacco et al.| (2008 1. 

Comparing the measurements of this work (see Paper II) and/or with those of 
^ Sources detected in our XMM-Newton survey (see Section 2.1). 



Sacco et al. 



(20081 



els depending on the resolution and the wavelength coverage of 
each campaign. As an example, for the medium resolution LRIS 
sample, we have only been able to study the bluest component 
of the optical K I doublet (the wavelength coverage of the instru- 
mental set-up did not reach the other component). We must note 
anyway that the EW of this component measured on the tem- 
plates has a very strong dependence on the luminosity class, and 
therefore we think that our estimated values listed in Table|6]are 
reliable. 

As mentioned before, we have also used both, templates and 
models (Allard et al.^2003) , to assess membership. In Fig|9]we 
show the cases of LOril35 and LOril46: with estimated tem- 
peratures of 3000 and 2800 K respectively, we have normalized 
both: the science spectra, and models for those temperatures and 
different values of \og(g). Even though the signal to noise ratio 
was not too high, it is very clear that these sources must have a 
value of \og{g) lower than 4.0 (and higher than giant-like values, 
where these doublets are barely detectable). 

In Fig. 10 we plot the effective temperature (in the range 



where the method is applicable, see Sec. 5.3 i of members and 
rejected candidates to CoUinder 69 against their measured Na I 
EW (the /18200 doublet). We have also included the measure- 
ments from Maxt ed et al.| ([2008 ) and highlighted in red those 
sources showing X-ray emission and/or sources that are classi- 
fied as active accretors in Paper II. LOril 15, a "peculiar" source 
(see further on, in the text), is marked as well. As we did with 
the Li I measurements, we have compared the values measured 
by us with those from Ma xted et al.| ( 2008| l for the sources in 
common. We find that ~35% of the common sources show Na I 
variability (differences not compatible within the error bars for 
both measurements). The majority of these sources are either 
Class II sources or X-ray emitters and therefore this variability 
must be mainly related either to activity or accretion processes. 

This diagram may suggest also that in the substellar domain 
(the boundary is located at TefF~3150 K for 5 Myr according to 
the SIESS + COND isochrone) active accretors show lower Na 
I EWs than non-accreting brown dwarfs. A first explanation for 
this trend would be that our measurements of Na I EWs are af- 
fected by veiling for the accreting brown dwarfs; this is not the 
case, since, as it is described in Paper II, even for the sources 
with the most intense Ha emission, the veiling is minimal for 
wavelengths redder than 7400 A. In any case, the sample is too 



small and the trend too shallow to allow us to perform any fur- 
ther analysis. Finally, the position of LOril 15 in the previously 
mentioned diagram seems suspicious (the Class II very low mass 
star clearly out of the general trend of the confirmed members), 
but its membership is discussed in the following section. 

5.3. Peculiar Sources: 

We distinguish in this section among four different groups: 
sources for which we have information about the Li absorption; 
objects for which we have information about other alkali lines 
(and/or emission lines and infrared excess); targets for which, 
due to the lack of high resolution spectra and its "hot" nature, 
we have had to use different criteria than alkali lines to study 
membership; and, finally, the coldest sample of sources only ob- 
served with SUBARU: 

1) Sources with Li information: 

LOri044, LOri046, LOri049, LOri052: AU these sources 
(except LOri052) have been classified as non-members based 
on the absence or very marginal detection of Li I. In the cases 
of LOri044, LOri046 , LOri 049, LOri052, either [Maxted et aL] 
(2008) orjSaccoetaL (2008 1 also classify them as non-members 
based on their radial velocities. In our FLAMES spectra of 
LOri052, a marginal detection of Li I in absorption with EW of 
~0.28 A was measured; this value is too low compared to other 
members with similar spectral type (even taking into account 
the general spread already commented on); therefore, together 
with the radial velocity determination of Maxted et al. ( 2008| l 
and our alkali lines study, we classify this source as probable 
non-member 

LOri007: No significant Ha emission was measured on the 
spectra of this source (observed with TWIN and with FLAMES), 
and although the equivalent width values of the Na I doublet is 
consistent with membership (see Table [6|l, no Li I in absorption 
has been detected (see Table[6]). Therefore we provide a classifi- 
cation of non-member 

DM065, DM070, DM061, LOriOOl, DM062, DM016, 
LOri026 (DM012), LOri038 (DM002): with s pectral types 
from K9 to M3: confirmed by D olan & Mathieu] ( |T999 ) via Li 
absorption. Besides, their positions in the HR diagram are com- 
patible with membership and most of them are X-ray emitters. 
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Fig. 9. Example of the surface gravity estimation for two NIRSPEC spectra by comparison with theoretical models. 
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Fig. 7. Lithium equivalent width (EW, in A) versus the spectral 
type for different sets of spectroscopically confirmed members 
of Collinder 69. Top: New members reported in this work where 
accretors (based on the Ha emission, see Paper 11 for details) are 
highlighted in red, and Class II sources (based on IRAC data, see 



Sigma Ori memb 
069 memb 




[Barrado y Navascues et al 



2007 and Paper II) and spectroscopic 



binaries (from |Sacco et al. 2008| and |Maxted et al.|2008| l are sur 



rounded by open circles and squares, respectively. Middle: New 
members showed as filled black circles, while members from the 



literature ( 


Dolan & Mathieu||1999 f2Ml 


jBarrado y Navascues 


et al.||2004 


c 2007|l are shown as smaller 


grey circles. Bottom: 



Comparison of Collinder 69 members (joining the samples from 
the middle panel), displayed as filled black dots, with Sigma 



Orionis low mass stars from |Zapatero Osorio et al.| (2002), rep- 
resented with star symbols. Regarding all the panels: The solid 
line traces the upper envelope of the values measured in older 
clusters such as IC2391, IC2602, the Pleiades and M35; and the 
short-dashed blue lines corresponds to the cosmic abundances - 
A(Li) = 3.1 - from gravities of log g = 4.5 and 4.0, respectively 
(curves of growth from |Zapatero Osorio et al.|2002| l. 



• Accrctor 

IR excess 

1 I Binary 




Sacco et al. 



(20081 or 



Fig. 8. Lithium equivalent width listed in 

|Dolan & Mathieul ( |1999| |2001| l versus the values measured in 
this work, in A. We highlight especial sources in red (accretors), 
or surrounding them with a circumference (objects showing in- 
frared excess) or a square (binary systems according to [Sacco 
et al. 2008 ; M axted et al.|2008| l. Vertical or horizontal sohd hues 
represent the range in variation when more than one measure- 
ment existent for one of the axes. 



LOrillS: Shows a deeper Na I absorption (LRIS spectrum) 
than expected for its spectral type (M5) but a Li EW (FLAMES 
spectrum) perfectly compatible with membership. Besides, the 
SED of LOril 15, exhibits a clear infrared excess at IRAC wave- 
length suggesting that the source is surrounded by an optically 
thick disk. It also shows quite strong and variable Ha emission 
(EW of -12.5 and -9 in our FLAMES and LRIS spectra), very 
close to the limit where pure activity cannot explain this emis- 
sion. All this clear signs of youth allow us to classify it as mem- 
ber and propose as possible explanation for the deep Na I ab- 
sorption some phenomenon related with activity. 
2) Other alkali lines and extra information: 
LOri036: shows faint Ha emission, but the strength of the 
alkali lines suggests that this object it older than the cluster mem- 
bers. Furthermore |Sacco et al.| ( [2"008| l found a RV value not fully 
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3000 
T (K) 



Fig. 10. Teff vs EW(Nal) for members and non-members of 
Collinder 69. Rejected candidates are displayed as crosses, con- 
firmed members as solid dots, where: red ones highlight sources 
classified as active accretors and/or detected on our X-rays sur- 
vey. Large open circles and squares surround sources showing 
infrared excess or classified as binaries in either Sacco et al.l 
( 2008| l or Maxted et al. ( 2008[ l . For those sources showing large 
variability in the measured Na I equivalent widths (measure- 
ments from [Maxted et al. [2008 [ and this work) we have included 
vertical bars connecting those measurements. Note the unusual 
position of LOril 15 (discussed in Section [53 i 




T _ (K) 



Fig. 11. HR diagram for the ~200 spectroscopically confirmed 
memb e rs by [Dolan & Mathieu[ ( [T9991 |200T] i; [Maxted et aT] 
( [2008^ ; [Sacco et alj (|2008| or this work. We hav e over-plotted 
several isochrones and evolutionary tracks from Baraffe et al. 
( [T998l l. Symbols as in Fig.[TO] 



compatible with membership. For these two reasons we classify 
this object as possible non-member. 

LOri046, LOri049: None of these sources shows He emis- 
sion; both have very low gravity according to the alkali absorp- 
tion strengths; and both had indeed been classified already as 
possible non-members by Sacco et al. (2008]) because of their 
measured radial velocity. We classify as no members. 

LOrillO, LOril33, LOril41, LOril47, LOrilSl, 
LOril54, LOril65: All these sources show surface gravities 
(based on the measured EW(Nal), see Table |6| too large 
compared with other cluster members. For the particular case 
of LOrillO, no Li I was detected on our LRJS spectra, and 
[Maxted et al.[ ( [20()8 ) also measured Na I EWs inconsistent with 
membership and confirmed that the radial velocities derived 
ruled out membership too. Regarding LOril33, LOril51 and 
LOril65, [Barrado y Navascues et aL ( 2004c| l already classified 
these objects as doubtful candidate members based on their 
photometric properties 

3) Sources with spectral type earlier than M3 and no 
Lithium information: 

As verified by [Slesnick et"aL] ( |2006| , the Na I -8200 
A doublet strength saturates for spectral types earlier than 
M2. Furthermore, the differences among dwarfs, giants and 
young sources become very subtle already at M3 spectral types. 
Therefore, for these "hotter" sources where we did not have 
spectra with resolution high enough to measure Li, we have used 
different sets of criteria to assess membership: 

- Non members: 

C69-IRAC-010 and C69-IRAC-008 with spectial types F9 
and K3, show Ha in absorption, no infrared excess at the IRAC 
wavelength range and a position in the HR diagram well below 
the 100 Myr isochrone, therefore we discard them as members. 

- Members: 

C69XE-009 (K2 spectral type) shows a very intense Ha 
emission and the estimated position in the HR diagram is per- 
fectly compatible with membership. 

C69-IRAC-001 and C69-IRAC-002 (Ml. 5 and M3 spectial 
types, respectively): Both are class II sources with thick disks 
according to the IRAC photometry and their positions in the HR 
diagram are compatible with membership. 

LOri024, LOri061, LOri048 and LOri062 (spectra l types 
from Ml. 5 to M3): All of them have been confirmed by [Sacco] 



etal. ( 2008 1 and/or Maxted et al. ( 2008 1 via radial velocity (and 



in all the cases, their positions in the HR diagram are compatible 
with membership). 

C69-IRAC-007 and C69-IRAC-005 (M2.5 and M3): Both 
sources are active accretors (see Paper II) with thick disks and 
positions in the HR diagram compatible with membership. 

C69 XE-072 (M3 sp ectral type): Detected in our X-ray sur- 
vey ( [Barrado et al.pon ) and with an estimated position in the 
HR diagram compatible with membership. 

4) Subaru L type sources: The resolution of our 
Subaru/IRCS spectra was too low to resolve any alkali line. 
However, given the confirmed spectral type, their luminosities 
are not compatible with objects being much closer or much fur- 
ther than Lambda Orionis and therefore we assess Y? member- 
ship only expecting to confirm definitely membership by means 
of, for example, a proper motion study in a longer time baseline 
than the one provided by the images from which these candi- 
dates were selected (we only have one set of Subaru images; see 
iMorales-Cal der6n|2008| l 

To summarize and conclude the section, in Table ]6] we have 
included a column with the final membership assigned to each 
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Fig. 12. Na 1 /18200 doublet detail for some of the particular 
sources (all but LOrillS discarded as members) with CAFOS 
spectra. For each case we have over-plotted in red a confirmed 
member with the same spectral type, and in blue a field dwarf 
(again with the same spectral type and obtained in the same cam- 
paign). Note that for spectral types M4.75 and M5.25 we have 
averaged spectra of field dwarfs with spectral types M4.5 and 
M5.0, M5.0 and M5.5, respectively. 



source. We consider members those labeled as Y?: objects show- 
ing some peculiarity already discussed but still considered as 
very good candidates for membership; and Y: spectroscopically 
confirmed members. In the latter case, we have included the 
reference of the work were the spectroscopic confirmation was 
achieved. 

There are ~240 members from which 225 have been labeled 
as Y? or Y (36 of them are confirmed as members of Collinder 
69 for the first time in this work). These confirmed members are 



displayed in an HR diagram in Fig 1 1 We must note that there 
is a large dispersion on the diagram for the hotter part of the 
sample (Teff larger than ~3900 K), where almost 100% of the 
candidates were confirmed by Dolan & Mathieu (1999 2001[ l. 
However, if we take a closer look to the M (and cooler) popula- 
tion of Collinder 69 (where this work is focused); the best fitting 
isochrone is that of 5 Myr with an upper limit value of 20 Myr 
recovering 95% of the confirmed members. The dispersion in 
this narrower area of the diagram, although less significant than 
the one affecting hotter sources, is mainly explained by objects 
harboring disks and has been previously addressed in |Bayo et al.| 
( |2008l l. 

Focusing only on the spectroscopically confirmed members 
(from this work and the Uterature), we find a brown dwarf to star 



ratio (/?ss as defined by Briceno et al. 2002) of 0.06. This value 
is similar to that found for Taurus by Briceno et al.|p002)l but 
significantly smaller than the revised one by ,Guieu et al.^2006| l 



or the value reported for the ONC (Kroupa & Bouvier 



2003] l. 



This fact can either pointing towards environmental effects in 
the substellar formation efficiency (as suggested previously in, 
for example Briceflo et al.| 2002 ; Kroupa & Bouvier 2003 i or 
towards pure limitations in the spatial extension of our surveys. 



6. Mass segregation: ejection mechanism. 

As mentioned in Sec.[T[ one of the main goals of our long term 
project on the LOSFR is to investigate which of the currently 
proposed mechanisms of formation of brown dwarfs agrees best 
with the observations. 



According to the ejection scenario by Reipurth & Clarke 
( |2001 ), brown dwarfs of an association such as Collider 69 
should be located mostly around the exterior parts of the clus- 
ter. If we assume an escape velocity ~1 km/s (as for the case 
of Taurus, Kroupa & Bouvierl|2003 1, after ~5 Myr, the ejected 
brown dwarfs should still have been recovered by some of our 
photometric surveys (a ~35' distance is well covered by our 
SUBARU field and marginally by the Spitzer/IRAC one). 

To test this hypothesis, we have studied the spatial distribu- 
tion of both stars and brown dwarfs in Collinder 69. In Fig. 13 
we plot two diagrams corresponding to the distribution of stel- 
lar and substellar confirmed and candidate members to Collinder 
69 (left panel for the whole set of confirmed members and good 
photometric candidates, and the right one only showing the spec- 
troscopically confirmed members). In both panels stellar sources 
are displayed as small black stars whilst substellar objects are 
highlighted with red dots. Even though what we can see in this 
figure is nothing but a projection, any structured grouping of the 
substellar population is far from obvious. 

From a more quantitative point of view, we tried to com- 
pare the distribution of stars and brown dwarfs with a two di- 
mensional KS test, but because of the low numbers we are deal- 
ing with, the results were inconclusive. Trying a more simplistic 
approach, since we are mainly interested in the radial distribu- 
tion of brown dwarfs, we built the cumulative radial distribution 
of stellar and substellar sources (good candidate members and 
spectroscopically confirmed members, as in Fig. 13 i. 



Although simple, this analysis had to be done carefully tak- 
ing into account again the low numbers we are dealing with. 
To obtain the cumulative distributions in a robust way, we per- 
formed ~60 different histograms using both multi-binning and 
multi-starting-point approaches and then smoothed the resulting 
function to make sure that the features found are bin and starting 
point independent. 



We show the results of this method in Fig. 14 where the 
data-points of the individual histograms are drawn as black dots, 
and the smoothed final function is highlighted in red. It seems 
quite clear that there is no obvious differences in the distri- 
bution of stellar members of Collinder 69 (supposedly formed 
according to the classical paradigm) and the substellar ones. 
Furthermore; if we compare the distribution of brown dwarfs 
with that of objects that are clearly stars (masses larger than 0.5 
M0), the results remain the same. 

The only evident feature in these radial distributions is a void 
of both stars and brown dwarfs at ~8-12 arcmin distance from 
A Ori (projected ~1-1.4 pc). We have highlighted this void in 
Fig. 13 were we see no correlation of this area with differences 
in dust densities as traced by the IRAS contours (the solid black 
lines). 

To try to quantify the physical relevance of this void, we have 
produced 100 synthetic 2D normal populations of sources. For 
each case we have produced three different sets of data: stars, 
brown dwarfs, and star + brown dwarfs. To draw each set, we 
have used the same number of objects that we have spectroscop- 
ically confirmed for each category in Collinder 69. The area from 
which the objects are drawn is the same area covered by our pho- 
tometric surveys. 
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Fig. 13. Spatial distribution of: left panel: good photometric candidates and right panel: spectroscopically confirmed members. O 
and B populations are highlighted with large blue (four pointed) stars, members with masses above the H burning limit are displayed 
as small black (five pointed) stars and brown dwarfs as filled red dots. In both figures it is clear that brown dwarfs are not located 
at the edges of the cluster (as the ejection mechanism would predict). We have shaded in grey the ring corresponding to the void in 



the radial profile displayed on Fig. 14 



We have run a one dimensional KS test on the observed ra- 
dial distribution compared to the synthetic ones (given the low 
number of brown dwarfs we have worked exclusively with the 
"only stars" and "stars + brown dwarfs" sets) and the result 
is that the probability of the observed radial distribution being 
drawn from a normal 2D population is negligible. Furthermore, 
only 7% of the synthetic radial distributions show voids similar 
in width (but closer to the edges of the cluster) to the one present 
in our observed distribution. Therefore we are confident that this 
void is not just a statistical artifact. 

Another interesting issue also not supporting the ejection 
scenario is the fact that Collinder 69 seems more extended in the 
stellar population than in the substellar one (the vertical solid 
line in Fig. [T4| marks the distance after which no more substel- 
lar members are found). This fact would indicate that the cluster 
has begun to be dynamically relaxed, so low-mass members are 
falling to the gravitational well. An alternative, given the very 
young age of the cluster, would be that the brown dwarfs might 
just have formed preferably closer to the center As a note regard- 
ing the possible effect of the spatial coverage of our surveys, we 
want to point out that even though our SUBARU field reaches 
further separations from Lambda Orionis than the other surveys, 
we do not find brown dwarf candidates more distant than ~20 
arcmin. 

In any case, we must remember once again that we are deal- 
ing with small numbers and that our photometric surveys are not 
wide enough to recover the possible ejected brown dwarfs if their 
velocity escape would be twice the one estimated for Taurus. 
Therefore, the most we can conclude from our analysis is that 
the ejection mechanism of formation does not seem to be able 
to explain the observed properties of Collinder 69 for velocity 
escapes of ~1 km/s or lower. 



7. The Initial Mass Function for spectroscopically 
confirmed members 

To derive the mass of the members of Collinder 69, we have used 
the estimated Tejf and Lboi from the SED fit along with a com- 



posite 5 Myr isochrone ( jSiess et al. 2000_ 4- COND from ,Baraffe 
|etal.|2003] l. The average mass resulting from these two estimates 
for each object is provided in Table|6] As can be seen in Fig. 15 
there is a quite large dispersion when comparing the estimates. 
If we separate our members according to the ratio between the 
two estimated masses into "large dispersion" (those where one 
of the masses is larger than 1.5 times the other one) and "accept- 
able" dispersion (the ratio between the largest and the lowest 
mass lower than 1.5); we find that ~3Q% of the so-called "large 
dispersion" sources have been classified as Class II according to 
their IRAC colours, while the same class in the acceptable dis- 
persion sample decreases to 15% (this is not a question of scale 
since both samples have a similar number of sources). In fact, 
~65% of the sources undergoing active accretion fall into the 
large dispersion sample. Therefore, our mass estimates are quite 
sensitive to the presence of disks (a caveat that does not surprise 
us, as has already been discussed in |Bayo et al.||2008| . To take 
into account this uncertainty, in this section we will present all 
the calculations for both estimates of the mass of the members 
and for the average masses. 

In the left panels of Fig. [16] we present the derived Initial 
Mass Functions (according to the different procedures followed 
to derive the masses of the sources) for Collinder 69. For this 
first case, we have considered only spectroscopically confirmed 
members confined within the intersection of the CFHT, Spitzer 
and XMM fields of view (see Fig. [TJ. For the assembly of each 
IMF, in the range of masses covered by this work, we have pro- 
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Fig. 14. Radial distribution of stars and brown dwarfs in 
Collinder 69: distance to A Ori vs. cumulative distribution of 
sources normalized to the radius where no more brown dwarfs 
are found in Collinder 69. As in the previous figure we include 
two sets of panels: on the left side candidates and confirmed 
members (top panel for stars and bottom for substellar objects) 
and on the right panel only spectroscopically confirmed mem- 
bers (as in the previous case, top for stars and bottom for substel- 
lar sources). Note the obvious void (grey shaded area) of brown 
dwarfs (and not so obvious but still present for stars too) from 
~8-12 arcmin distance to A Ori. 



ceeded in a similar manner as in the previous section with the 
multi-binning, multi-starting-point approach. Also, to enlarge 
the mass range coverage, we have included data from Murdin 



& Pension ( 1977 1 (scaling to the respective areas) as open blue 



circles. The last point at the massive end correspond to a possi- 



ble SN (open blue triangle); see discussion on Dolan & Mathieu 
([1999). 

We must note that the previously described area (the inter- 
section of the three FoVs), leaves the central part of the cluster 
(where the 08 111 binary star, A Ori, is located) aside. In order 
to study the possible effects of this fact, in the right panels of 
Fig. 16 we show another set of IMFs (again three for the differ- 
ent estimated masses), this time focused on the CFHT field of 
view (scaling the XMM members to the corresponding area). 

Regardless of the parameter used to infer the mass of our 
confirmed and candidate members, for masses above 0.65 Mq, 
the power law index of the IMF is similar to Salpeter's value, 
being much smaller for lower masses. In particular, the fitted 
slopes in the mass range 25-0.70 Mq vary from 1.73 to 1.89; 
and those for the low and very low mass range, 0.70-0.01 Mq, 
correspond to much lower values; from 0.16 to 0.38 (see Table|5] 
for a summary of the slopes). 

As we have mentioned before, the two areas for which we 
have derived the IMFs (the CFHT survey coverage and the inter- 
section of CFHT, Spitzer and XMM) comprise populations that 
show distinct trends in terms of spatial distribution (the XMM 
coverage leaves a side a relatively large population of Class II 
sources). Surprisingly, this fact does not seem to affect the shape 
of the IMF since the slopes shown in Fig. 16 agree among them- 
selves. 

Our estimated slopes are systematically lower (but in a simi- 
lar range) than those derived for photometrically compiled IMFs 




® 
I 



0.10 



Fig. 15. Relationship between the masses derived from Tgff and 
Lboi. Members (spectroscopically confirmed) and candidates 
(photometric) are displayed as solid dots. As in previous figures 
we have highlighted peculiar sources (Class II, binaries and ac- 
cretors). 



Table 5. Fitted slopes for the IMF of Collinder 69 taking into 
account the different methods to estimate the mass of the con- 
firmed members (see text for further details). We also include 
the slopes estimated for low mass stars (the mass ranges are not 
always the same ones but roughly M* < 0.8 Mq) in other works 
for different open clusters. 



Estimated a for 
the mass range: 


25-0.65 Mo 




0.65-0.01 Mo 


C69 IMF (I) 


1.73 




0.18 


C69 IMF (II) 


1.81 




0.38 


C69 IMF (III) 


1.88 




0.31 


C69 IMF (IV) 


1.79 




0.16 


C69 IMF (V) 


1.85 




0.38 


C69 IMF (VI) 


1.85 




0.33 


SOri' 




0.8 




Apgj.2.3 




0.6 




M35* 




!l)-(-0.8J 


5) 


Pleiades^ 




0.6 





S lope estimated wit h data from: 
'| B6iaretal.1j200l|l 



^ Barrado y Navascues et al. 1 2001a 
^ Barrado y Navascues et al. (zOOST 
Barrado y Navascues e t al.^^2001b) l 
^ |Bouvieretal.| ( |1998^ ' 



for several open clusters: SOri, APer, M35, and the Pleiades (see 
Table [5]for the exact numbers and references). For the low mass 
stars/substellar domain, the differences are not too large, regard- 
less of the total mass, the environments and the age of each as- 
sociation. The only exception to this is the case of M35 -which 
corresponds to the core of the cluster-, a very rich association. 
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about 150 Myr old, where important dynamical evolution might 
have taken place. 



Fiducial lognorm MF: Mc=0.2a, 0=0.52 (Pleiades, Moraux et al. 2003) 
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Fig. 16. IMFs derived for Collinder 69 taking into account three 
different determination of the mass for each members and two 
different areas within the cluster. Left panels: Area confined 
within the intersection of XMM, CFHT and IRAC fields of view. 
Top: IMF (I). Masses derived combining the estimated Tejj and 
the composite 5 Myr isochrone (SIESS + COND). Middle: IMF 
(II). Masses derived combining the estimated Lboi and the com- 
posite 5 Myr isochrone (SIESS + COND). Bottom: IMF (III). 
Masses derived averaging the masses from the previous two fig- 
ures. Right panels: Area confined within the CFHT field of view 
(including the central part of the cluster). Top: IMF (IV). Same 
method to estimate masses as in IMF (I). Middle: IMF(V). Same 
method to estimate masses as in IMF (II). Bottom: IMF (VI). 
Same method to estimate masses as in IMF (III). 



In particular, if we only compare with similar age clusters, 
as we do in Fig. 17 we can see that our IMF (in this case we 
have used the mean mass determination for the intersection of 
the three fields of view) looks very similar to that of, for exam- 
ple, NGC6611 with, in our case, a cut-off at ~0.016Mq, were 
we think (because of the completeness of our photometric sur- 
veys) that we have reached the minimum mass for members of 
CoUinder 69. 



8. Conclusions 

We have analyzed ~170 optical and near infrared spectra with a 
wide range of resolutions of candidate members to the Collinder 
69 cluster Based on different criteria regarding molecular ab- 
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Fig. 17. Log-log representation of the IMF for different young 
associations. We compare our very complete spectroscopically 
confirmed IMF (in pink) with photometric ones derived for as- 
sociations with simi lar ages (~2-5 Myr, |Luhman et al.| 2007} 



Pliveira et al.|2009t|Caballero et al.|2009| |Lodieu et al.|2007| l. 

The compilation of the data corresponding to clusters other than 
Collinder 69 are courtesy of J. Bouvier and E. Moraux. 



sorption bands, we have provided spectral types for all the 
sources. Using alkali lines as youth indicators, we have con- 
firmed 90 members (and 9 possible members). For those sources 
of our survey overlapping with radial velocity surveys (by Sacco 
et al.|200^ Maxted et al.||2008] l, we obtain very similar results 
regarding membership. A summary of the sources analyzed in 
this work can be found in Table |6] 

We have confirmed the cool nature of the lowest mass can- 
didate members to Collinder 69 reported so far With derived 
spectral types L0-L2 corresponding to effective temperatures of 
~ 2000 K (according to the spectroscopic temperature scales 
derived by Basri et al. 2000), if members, they have estimated 
masses, following the 5 Myr DUSTY isochrone by Chabrier 
et al.|2000i as low as 0.016 Mq (-16.5 Mjupitei). 



iLpi 



e have built the most complete IMF based on spectroscop- 
ically confirmed members for a young association and when 
compared with others (photometric IMFs), no significant differ- 
ences in the low mass and substellar domain are identified. 

With our study of the spatial distribution of confirmed stel- 
lar and substellar members of Collinder 69, we provide some 
caveats to the ejection formation scenario for brown dwarfs, in 
order for it to be feasible in this region. 
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Appendix A: Automatic line characterization 
procedure and resolution effects. 

In short; one of the most difficult steps to automate when mea- 
suring Hne properties is the determination of the continuum on 
top of (or below) which the line lies. Even when this task is 
performed manually, different astronomers might select different 
points as the bases of the continuum when, for example, facing 
a noisy spectrum or a line embedded in a molecular band. The 
code that we have developed tries to solve this problem proceed- 
ing in an iterative and consistent way. 

First of all it locates where the actual peak of the line is; in 
this step the wavelength calibration accuracy is taken into ac- 
count providing limits in the wavelength coverage where the 
search for the peak of the line is carried out. Once this task is 
accomplished, two regions are defined (to the right and to the 
left of the estimated central wavelength) the width of the regions 
is fixed to 100 points to be statistically significant (and therefore 
the width in wavelength units would depend on the resolution of 
the spectra). A linear fit using those regions is performed (red 
dotted line in Fig A.l i. With this attempt to define a continiuum 



6562,80 6562.52 



(the linear fit), the code calculates a first guess of the FWHM 



(red cross in Fig A. 1 



For the second iteration, the code consider s that the l ine is 
a gaussian, and, therefore, FWHM - 2 x -y/2 x log(2) x cr. 
Assuming that at a distance of +10cr one should be outside the 
line, two new regions on the spectra are selected: starting at 
ilOcr from the wavelength of the peak and ending at the same 
limits as before. These two regions are close to the edges of the 
line unless the line has a very wide double peaked structure. 

At this point the code will perform two linear fits in a sequen- 
tial way; a second continuum is derived with a linear fit to these 
new regions of the spectra close to the edges of the line, and 
the third one with another lineal fitting, but this time consider- 
ing only those data-points that differ from the second continuum 
less than once the dispersion of the difference between the actual 



spectra and this second continuum (the blue dots in Fig A. 1 rep 



resent those data-points considered to define a third continuum 
and the blue dashed line is the resulting one). 

Once the third guess on the locus of the continuum is es- 
timated, the code calculates a second iteration of the FWHM 
(blue cross in Fig A.l i. Assuming one more time that the line 



can be described with a gaussian (and calculating the associated 
cr), it considers three pairs of lambdas in the spectra to make the 
final measurements (pairs located at ±3cr, +4cr and ±5cr from 
the lambda of the peak). As a final refinement, for each selected 
pair of wavelengths; for example Ai and A2 corresponding to 
'tpeak ± 3cr the code checks whether the linear continuum defined 
by {A\,¥{A\)), (/l2,F(/l2)) intersects the wings of the line; when 
this is the case (for example the upper-left panel of Fig A. 1 1, the 
code will use the third continuum to perform the measurements. 

Prior to provide these measurements, the code subtracts the 
instrumental profile from the values estimated for the FWHM 
and FWi()% considering that the FWHM of the convolution 

(Gmeasured) of tWO gaUSSianS (Ginstr, Gline) foUoWS: 

FWHM(G_ed) = 7FWHM(Gi„,,,)2 + FWHM(Gu„e)' (A.l) 

and that FWHM and FWi()% relate accordig to: FWio% - 
ViSxFWHM 

The final product is the mean of the three measurements for 
each parameter (FWHM, FWio% and EW) and their correspond- 
ing standard deviation cr. 




LOriOJ' LI 6707.80 6707.52 
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Fig. A.l. Example of the output of the code designed to mea- 
sure EW, FWHM and FWio% of the different lines detected on 
the spectra in an automatic manner. In this case we show four of 
the measured lines on the averaged MIKE spectra obtained for 
LOriOSl: Ha at 6562.80A, He I at 6678. 15A, Li I at 6707.8Aand 
S II at 67 17. OA. For each case the canonical wavelength for the 
measured line and the actual wavelength where the line has been 
found in our spectra are displayed o the title of the plot (all of 
them within the errors on wavelength calibration). The iterative 
process in the continiuum determination is displayed using dif- 
ferent colours: red for the first iteration (including the first guess 
for the half maximum flux value, red cross), blue for the sec- 
ond one (filled blue dots on the data-points used to refine the 
first continiuum, blue cross showing the second half maximum 
estimation) and the three green lines (filled, dotted and dashed) 
representing the final determined continuum. 



Once we had an automatic procedure to characterize the lines 
present in our spectra, we had to address the effect that the wide 
dispersion in resolution might introduce in our measurements. 
We proceeded in the following manner: We degraded one of our 
FLAMES spectra (from LOri038, a M3, accreting. Class II, con- 
firmed member) to different resolutions (including those listed in 
Table [T}. We used then (on those degraded spectra) our code to 
measure FHWM, EW and FWio% for th e Hce (emission) and Li 



A.2 



we provide those values 



/16708 A lines (absorption). In Fig. 
(EW and FWio% for Wa and EW for Li I 216708 A) as a func- 
tion of the resolution of the spectra. The length of the y-axes has 
been fixed by purpose for each plot: in the first panel (from left 
to right), it represents the variation of the instrumental response 
(the FHWM measured on the respective arc adapted for each 
value of the resolution) among the resolutions considered. In the 
cases of the middle and right-side panels, the length of this axes 
provides an idea of the variation within members of Collinder 
69 for the corresponding equivalent width. Note that in the case 
of Hot there is clearly no significant dependence on the measure- 
ments made at different resolutions (the variations lie within the 
error bars). For the case of the lithium, due to the intrinsic weak- 
ness of the line, the accuracy of the measurement (meaning the 
error-bars) lowers with the spectral resolution, but still it seems 
perfectly reasonable to compare measurements at R ~ 2000 and 
R ~ 8000. In the same figure we have highlighted the range of 
resolutions where Li is not detectable anymore. A detailed view 
of the "degeneration" of the lines for this example is provided in 
fig.|A3] 
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Fig. A.2. Relationship between different measurements of the Hor and Li I hne profiles and the resolution os the spectra where the 
measurement has been performed. 
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Fig. A.3. Sequence composed by a detailed view of the evolution of the profiles of the Ha (emission) and Li I (absorption) lines with 
the variation of the spectral resolution. The first panel (from left to right) corresponds to the original resolution of our FLAMES 
spectrum of LOri038, a M3 member of C69 (R ~8000). In the last panel we have degraded the FLAMES spectra to a resolution 
similar to that of our low-resolution campaign (R ~800 obtained with the B&C spectrograph). As can be seen in this figure and in 
Fig. A. 2 the lithium's equivalent width can be measured down to resolutions of the order of 1250 (for a good signal to noise ratio as 
it is the case of our spectrum of LOri038). The linestyles and color code are those explained on the text. 
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